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Abstract. We consider travelling periodic and quasi-periodic wave 
solutions of a set of coupled nonlinear Schrodinger equations. In fi- 
bre optics these equations can be used to model single mode fibers 
with strong birefingence and two- mode optical fibers. Recently 
these equations appear as model, which describe pulse-pulse inter- 
action in wavelength-division-multiplexed channels of optical fiber 
transmission systems. Two phase quasi-periodic solutions for in- 
tegrable Manakov system are given in terms of two dimensional 
Kleinian functions. The reduction of quasi-periodic solutions to 
elliptic functions is discussed. New solutions in terms generalized 
Hermite polynomials, which are associated with two-gap Treibich- 
Verider potentials are found. 



1. Introduction 

We consider the system of two coupled nonlinear Schrodinger equa- 
tions 

iU t + U xx + (kUU* + xW*)U = 0, 

(1.1) 

iV t + v xx + ( x uw + pVV*)V = 0, 

where k, x, p are some constants. The integrability of this system is 



proven by Manakov |31[ only for the case k = x — P-, which we shall 
refer as Manakov system. 

The equations fll.lQ are important for a number of physical applica- 
tions when x is positive and all remaining constants equals to 1. For 



example for two-mode optical fibers x = 2 |T3[ and for propagation of 
two modes in fibers with strong birefringence x = § an d in the 
general case | < A < 2 for elliptical eigenmodes. The special value 
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X = 1 (Manakov system) corresponds to at least two possible cases, 
namelly the case of a purely electrostrictive nonlinearity or, in the el- 
liptical birefringence, when angle between the major and minor axes 
of the birefmgence ellipse is approximately 35°. The experimental ob- 
servation of Manakov solitons in crystals is reported in . Recently 
Manakov model appear in the Kerr-type approximation of photore- 
fractive crystals |HJ. The pulse-pulse collision between wavelength- 
division-multiplexed channels of optical fiber transmission systems are 
described with equations (Q x = % E§ S HI • 

General quasi-periodic solutions in terms of n-phase theta functions 
for integrable Manakov system are derived in |TJ , while a series of spe- 
cial solutions are given in || |38|, |36|, 37]. The authors of this paper 
discussed already quasi-periodic and periodic solutions associated to 
Lame and Treibich-Verdier potentials for nonintegrable system of cou- 
pled nonlinear Schrodinger equations in frames of a special ansatz [12]. 
We also mention the method of constructing elliptic finite-gap solutions 
of the stationary KdV and AKNS hierarchy, based on a theorem due 



to Picard, is proposed in |19L 20, UM and the method developed by 



Smirnov in series of publications, the review paper 1 41] and H3L H2 



In the present paper we investigate integrable Manakov system being 
restricted to the system integrable in terms of ultraelliptic functions by 
introducing special ansats, which was recently apllied by Porubov and 
Parker |37[ to analyse special classes of elliptic solutions of the Manakov 
system (k = x — P — 1)- More precisely, we seek solution of ( |1.1|) in 
the form 



1-2) U(x, t) = qi(x) exp < %a\t + iC\ 



dx 
q\{x) 



V(x, t) = q 2 (x) exp < ia 2 t + iC 2 J 



X 

dx 



ql{x) 1 ' 



where the functions q\, 2 {x) are supposed to be real and a,\,a 2 ,C\,C 2 
are real constants. Substituting (|1.2|) into ( |1 . 1| ) we reduce the system 
to the equations 



;i.3) lh? +Pq * + Xq ^ 2 ~ aiQl ~ ~qj = ° 



^ 2 ?2 I 3 i 2 ^2 



COUPLED NONLINEAR SCHRODINGER EQUATIONS 



3 



The system (|1.3|) is the natural hamiltonian two-particle system with 



the hamiltonian of the form 

, . % 1 9 1 o 1 C? 1 Co 

(L4) - 2^" 2^ + 2^ + 2^' 

where = dqi(t)/dt. 

These equations describe the motion of particles interacting with 
the quartic potential Aqf + Bq\q\ + Cq\ and perturbed by inverse 
squared potential. Nowdays four nontrivial cases of complete inte- 
grabilty are known for nonperturbed potential (i) A:B:C= 1:2:1, (ii) 
A:B:C= 1:12:16, (hi) A:B:C= 1:6:1, (iv) A:B:C= 1:6:8. Cases (i), (ii) 
and (hi) are separable in respectively ellipsoidal, paraboloidal and Carte- 
sian coordinates, while the case (iv) is separable in general sence ||39|| . 
The cases (ii) appears as one of the entries to polynomial hierarchy 
discussed in JL5| the cases (hi) and (iv) are proved to be canonically 
equivalent under the action of Miura map restricted to the stationary 
of coupled KdV systems associated with fourth order Lax operator 0. 
Moreover all the cases (i)-(iv) permit the deformation of the poten- 
tial by linear conbination of inverse squares and squares with certain 



limitations on the coefficients |T5|, y] . There are also known Lax repre- 
sentations for all these cases which yield hyperelliptic algebraic curves 
in the cases (i) and (ii) and 4-gonal curve in the cases (hi) and (iv). 

Although each from the system enumerated yield nontrivial classes 
of solutions of the system ( |Llp we shall discuss further only the case (i) . 
The integrability of this case and separability in ellipsoidal coordinates 
was proved by Wojciechowski (see also (PJ). We employ this 



result to integrate the system in terms of ultraelliptic functions (hyper- 
elliptic functions of the genus two curve) and then execute reduction 
of hyperelliptic functions to elliptic ones by imposing additional con- 
strains on the parameters of the system. 

The paper is organised as follows. In the first section we construct 
the Lax representation of the system, develop the genus two algebraic 
curve, which is associated to the system and reduce the problem to 
solution of the Jacobi inversion problem associated with genus two al- 
gebraic curve. In the section two develop the integration of the system 
in terms of Kleinian hyperelliptic functions which represent a natural 
generalization of Weierstrass elliptic functions to hyperelliptic curved 
of hidher genera; recently this realization of abelian functions was dis- 
cussed in H, |1(], [16|]. We explain in the section the outline of the 
Kleinian realization of hyperelliptic functions and give the principle 
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formulae for the case of genus two curve. In Section 4 we develop re- 
duction of Kleinian hyperelliptic function to elliptic functions in terms 
of Darboux coordinates for the curve admiting additional involution. 
In this way a quasiperodic solution in terms of elliptic functions is ob- 
tained. In the last section we construct a set of elliptic periodic solu- 
tions on the basis of application of spectral theory for the Hill equation 
with elliptic potential. 

2. Lax representation 

The system l:2:l(/t = x = p= l)isa completely integrable 
hamiltonian system 

-q-2 + {Qi + ? 2 )?i - aiQi ~ -J = 0, 

(2.1) 



+ (Qi + 9 2 )?2 - a 2 q 2 - 



with the Hamiltonian 



1C| 



(2.2) H = \ grf + \ (,? 4- - \a iq \ - '-a,,! + i| , 2 % , 



where the variables 92,^2) are the canonicaly conjugated vari- 

ables with respect to the standard Poisson bracket, {• ; ■}. 

This system permits the Lax representation (special case of Lax rep- 



resentation given in |28fl) 
dL(X) 



<k [M(A),L(A)], 

" V(X) U(X) \ _ ( I 



(2 - 3) L(A) " \W{\) -V(X) J' M ~\Q{\) 
is equvalent to the ( |2.1|) , where U(X), W(X), Q(X) have the form 



U(X) = _„(,)(! 

W(X) = a( A)(-A + f + f + i(p? + f) 
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1 / 2 C7 2 \ 1 
+ o(P2 + -T 



2 V li J ^ - a 2 
Q(X) = X-ql-ql 

where a(A) = (A — ai)(A — a 2 ). 

The Lax representation yields hyperelliptic curve V = (z/, A), 

det(L(A) - ~ul 2 ) = 0, 
where I2 be 2 x 2 unit matrix and is given explicitly as 



v 2 = 4{X-a 1 ){X-a 2 ){X 3 -X 2 {a 1 + a 2 ) + X{a 1 a 2 -H)-F) 
(2.4) - C?(A - a 2 ) 2 - C7 2 (A - ai ) 2 , 



where if is the hamiltonian ( |2.2j ) and the second independent integral 
of motion F, {F; H} = is given as 

1 1 11 

F = -^(Piq2-P2qi) 2 + ^(ql + q2)(aia 2 - -a 2 qf- -a^) 

1 2 1 2 1 (2a 2 - g|)Cg 1 (2 Ql - qf)Cl 

(2.5) - - Pl a 2 - - g 4 g • 

We remark, that the parameters Ci are linked with coordinates of the 
points (dj, z/(aj)) by the formula 

(2.6) C7 2 = --^£- i, j = l,2. 

a,- - a,- r 



Let us write the curve ( |2.4|) in the form 

(2.7) z/ 2 = 4A 5 + a 4 A 4 + a 3 A 3 + a 2 A 2 + a x X + a , 

where the moduli of the curve a, are expressible in terms of physical 
parameters - level of energy H and constants ai, a 2 , Ci, C7 2 as follows 

«4 = — 8(ai + a 2 ), 

a 3 = -4if + 4(a x + a 2 ) 2 + 8aia 2 , 

a 2 = 4H(ai + a 2 ) — AF — C 2 — C\ — 8aia 2 (ai + a 2 ), 

ai = AF(a x + a 2 ) - Aa x a 2 H + 2Cla 2 + 2Cla x + Aa\a 2 2 , 

ao = — Aa\a 2 F — C\a\ — C\a\. 

Let us define new coordinates /ii,/i 2 as zeros of the entry t/(A) to 
the Lax operator. Then 

/ 9S n 2 K - MlKgl - /^2) 2 (Q2-Ml)(Q2-/^2) 

(2.») g : = 2 , g 2 = I . 

a\ — a 2 a 2 — a\ 
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The definition of fix, fi2 m the combination with the Lax representation 
comes to the equations 

(2.9) Vi = V(jii) = -^U(jii), i = 1,2, 

wich can be transformed to the equations of the the formQ 

/■mi rw 

(2.10) ux = I dux+ dux, 

J a\ J 0,2 

/•Mi rfJ-2 

(2.11) u 2 = du 2 + du 



'2 



(12 



where dux } 2 denote independent canonical holomorphic differentials 
(2.12) dux = ~, du 2 = —. 

V V 

and u\ = a, u 2 = 2x + b with the constants a, b defining by the initial 
conditions. The integration of the problem is then reduces to the solv- 
ing of the Jacobi inversion problem associated with the curve, which 
consist in the expession of the symmetric functions of (fix, fi2, vi, v 2 ) as 
function of two complex variables (ux, u 2 ). 



3. Exact solutions in terms of Kleinian hyperelliptic 

functions 

In this section we give the trajectories of the system under consid- 
eration in terms of Kleinan hyperelliptic functions (see, e.g. [§], fL0|D, 
being associated with the algebraic curve of genus two (|2.7f) which can 
be also written in the form 

(3-1) v 2 =4nto(A-A,), 

where A« 7^ \ are branching points. At all real branching points the 
closed intervals [A 2 i_i, A 2 j], i = 0, ... 4 will be referred further as lacunae 
p9| |3^| . Let us equip the curve with a homology basis (01, a 2 ; bx, b 2 ) e 
Hi(y,1i) and fix the basis in the space of holomorphic differentials as 
in (|2.12|) . The associated canonical meromorphic differentials of the 
second kind dr T = (dri, dr 2 ) have the form 

(3.2, dr, = ^ + + dj , 2 = A^,, 

Av v 



x \w what follows we shall denote the integral bounds by the second coordinate 
of the curve V = V(v, A). Q 
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The 2x2 matrices of their periods, 

2u = I <j> dui J , 2a/ — I (p dui . 

\Ja k / k,l=l,2 \Jh J k,l=l,2 

2r) = dr ( , 2r( = ( (b dr t 

\Jak / k,l=l,2 \Jh / k,l=l,2 

satisfy the equations, 

uj uj — ujuj =0, i] to — i]to — — — 1 2 , r] r] —rjr] =0, 

which generalizes the Legendre relations between complete elliptic in- 
tegrals to the case g = 2. 

The fundamental a function in this case is a natural generalization 
of the Weierstrass elliptic a function and is defined as follows 

7T e 

a(u) = 



det ( 2a; ) V 1 1, . 



i<j<5 ( v" , » 

x exp {u T r/(2cj)" 1 ii} ^[e]((2cj)" 1 w|c; / cj" 1 ), 
where e 8 = 1 and #[£](i;|t) is the 9 function with an odd characteristic 
which is the characteristic of the vector of Riemann 



£i ^2 



e' e' 



-1 c 2 

constants, 

9[e](v\r) = ^ ex P 171 {( m + £) T T(m + e) + 2(« + e') T r(m + e)} . 

mez 2 

Alternatively the a function can be defined by its expansion near 
u = 

(3.3) <r(v) = Mi + — a 2 u\ - -w 2 + o(u 5 ) 

and the further terms can be computed with the help of bilinear dif- 
ferential equation ||. 

The a-function posses the following periodicity property: put 

E(m, m!) = rjm + rj'm', and f2(m, m') = cum + u>'m', 

where m, m' G Z n , then 

a[e\(z + 2fi(ra,ra'),u;,a/) = exp{2^ T (m, m!) (z + f2(m,ra'))} 

x exp{— mm T m' — 2me T m'}a[e](z, uj, uj') 

As modular function the Kleinian cx-function is invariant under the 
transformation of the symplectic group, what represents the important 
characteristic feature. 
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We introduce the Kleinian hyperelliptic functions as second logarith- 
mic derivatives 

d 2 d 2 
Pu(u) = -^lno-(w), p 12 (u) = - — du In <?{u), 

d 2 

p 22 (u) = -—]n<T{u). 

The multi-index symbols pijk e ^ c - are defined as logarithmic deriva- 
tives by the variable Ui,Uj,Uk on the corresponding indices i,j, k etc. 

The principal result of the theory is the formula of Klein, which reads 
in the case of genus two as follows 



- ' f /Li p/J.2 



fej=l 



Pki ( / du - / du - / du fi k 1 f i l 1 



(3.4) = y i = l,2, 
4(/i - Hi) 2 



where 

2 

(3.5) F(/Zi, /i 2 ) = ^/i2[ 2 «2r + «2r+l(>l + Afc)]- 

r=0 

By expanding these equalities in the vicinity of the infinity we obtain 
the complete set of the relations for hyperelliptic functions. 

The first group of the relations represents the solution of the Jacobi 
inversion problem in the form 

(3.6) A 2 - p22{u)X - pn{u) = 0, 

that is, the pair (/ii,/i 2 ) is the pair of roots of (|3.6|) . So we have 

(3.7) p22{u) = A*i + {J, 2 , Pn(u) = -/Ui/i 2 . 
The corresponding z/j is expressed as 

(3.8) Ui = p222{u)^i + pi22(w), i = l,2. 

The functions p22, pn are called basis functions. The function pn(u) 
can be also espressed as symmetric function of fii, p.2 and Vi,V2'- 

(3 ' 9) P11{U) = 4^-^ ' 

where F(/ii,/i 2 ) is given in (|3.5|) . 
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Further from (13.81) we have 



/ x ^1 — ^2 / x ~ 

P222{U) = , p 221 (U) 



P2U{U) 

(3.10) p ul {u) 
where 



fll — /i2 ^1 — ^2 

[x\v 2 - fJ%V X 



wMjH, H2) ~ V\^{\i2i Vi) 
4(/i!-/i 2 ) 3 



H2) = 4a + «i(3/ii + /i 2 ) + 2a 2 /ii(/Ui + /i 2 ) 

+ a^nlini + 3/i 2 ) + 4a 4 Aii/i 2 + 4/^/x 2 (3/xi + /i 2 ). 

The next group of the relations, which can be drived by the ex- 
panding of the equations (|3.4|) are the pairwise products of the pijk 
functions being expressed in terms of p 22 , pi 2 , pu and constants a s of 
the defining equation ( |3.1| ). We give here only basis equations 

PI22 = 4 P22 + 4pi 2 p 22 + « 4 p 22 + 4pn + a 3 p 22 + a 2 , 
P222P122 = 4pi 2 p 22 + 2p\ 2 - 2p u p 22 + a 4 p 12 p 22 

+ \&Zpl2 + 

P 2 i22 = 4p 22 pl 2 - 4p n p 12 + a 4 pl 2 - a . 

All such expressions may be rewritten in the form of an extended cubic 
relation as follows. For arbitrary l,k G C 4 the following formula is 
valid H 

(3.11) l T nn T k = --det ( ? T [ 

4 y k 

where tt t = (p 222 , —p 2 2i, P211, —pin) and H is the 4 x 4 matrix 



(3.12) H 



\ai a 2 + Ap u \a z + 2p 12 -2p 22 
-2pu i«3 + 2p 12 a 4 + Ap 22 2 
\ ~2p 12 ^ -2p 22 2 / 



The vector tt satisfies the equation Hit = 0, and so the functions 
P22, P12 and pu are related by the equation 

(3.13) det# = 0. 



The equation (|3.13|) defines the quartic Kummer surface K in C 3 f23" 



The next group of the equations, which is derived as the result of ex- 
pansion of the equalities ( |3.4| ) are the expressions of four index symbols 
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p ijH as quadrics in p {j 

(3.14) P2222 = 6^22 + | a 3 + "4^22 + 4pi2, 

(3.15) P2221 = 6p22pl2 + "4^12 - 2pn, 

(3.16) P2211 = 2p22pn + 4p?2 + §«3pi2- 

(3.17) p21U = 6pi 2 pll + « 2 pl2 - ^«ip22 - «0, 

(3.18) p ull = §p\ x - 3a p22 + «ipi2 + "2^11 - |«ott4 + |«i«3- 

These equations can be identified with completely integrable par- 
tial differential equations and dynamical systems, which are solved in 
terms of Abelian functions of hyperelliptic curve of genus two. In par- 
ticular, the first two equations represent the KdV hierarhy with "times" 

(*1,* 2 ) = (W2,«l) = OM), 

(3.19) X k+1 [\J]=KX k [U] 

where 1Z = d% — U + c — lU^d -1 , c = 04/12 is the Lenard recursion 
operator. The first two equations from the hierarchy are 

(3.20) U tl = U x , U i2 = ^(U xxx - 6U X U), 



the second equation is the KdV equation, which is obtained from Q3.14 ) 



as the result of differentiation by x = U2 and setting U = 2p 2 2 + 
0(4/6. The equation (|3.14j) plays role of the stationary equation in 
the hierarchy and is obtained as the result of action of the recursion 
operator. 

Let us introduce finaly the Baker- Akhiezer function, which in the 
frames of the formalizm developed is expressible in terms of the Kleinian 
cr-function as follows 



a ( £ du - u ] r /-a 



a it 



(3.21) \fr(A, u) = v "~ , ^exp<j / dr T n 



where A is arbitrary and u is the Abel image of arbitrary point (u\, fi\) x 
(^2, ^2) £ V x V . It is straighforward to show by the direct calcula- 
tion, being bases on the usage of the relations for three and four-index 
Kleinian p-functions that \I/(A, u) satisfy to the Schrodinger equation 



(3.22) - 2p 22 (u))*(A, u) = (\ + ^o 4 ) * (A, u) 

for all (is, /i). 

Now we are in position to write the solution of the of the system 
in terms of Kleinian cx-functions and identify the constants in terms of 
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the moduli of the curve. Using (fT7| ) , (|2.8| ) the solutions of ( |2.1|) have 



the following form in terms of Kleinian functions ^22(^)5 pi2(u) 
2 n a i ~ P22(u)a! - pn(u) 

di — CLi 

(3.23) ql = 2 ^-P22(")°2-M«) > 

where the vector u T = (a, 2x + b). 

4. Periodic solutions expressed in terms of elliptic 
functions of different moduli 

We consider in this section the reduction Jacobi (see e.q.[^] ) of 
hyperelliptic integrals to elliptic ones, when the hyperelliptic curve V 
has the form 

(4.1) w 2 = z(z - l)(z - a)(z - (3)(z - a/3) 



The curve (4.1) covers two-sheetedly two tori 

tt ± : V = (w, z) -> E± = (rj±, £±) 

(4.2) ^ = e ± (i-e±)(i-4e±) 

with Jacobi moduli 



(4.3) k 



± (l-a)(l-/3)' 



The covers 7r± are described by the formulae 

(4.4) , ± = -V(l-a)(l-/n ( ,_^ a g^ . 

(4.5) £ = = < 1 -°»' 1 -^ . 

(2 — a) (2 — pj 

The following formula is valid for the reduction of holomorphic hyper- 
elliptic differential to the elliptic ones: 

(4.6) ^ = -^(l- a )(l-(3)(z =F v 7 ^)^ 



77± w 

Suppose that the spectral curve ( |2.7| ) admits the symmetry of the 



(|4~1"| ) and apply the discussed reduction case to the problem. Then the 
equations of the Jacobi inversion problem ( 2.11|) can be rewritten in 
the form 
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2ua 



2u. 



(4.7) v^KME / {z-^)^- 

i=l Jz o W 

2 fZi A~ 

(4.8) ^(1 -/?)(!- «) E / (^+V / ^)- 

i=l ■'so W 

with (^,/ii) = (2tWj,Zj) and 

(4.9) w± = - y/(l-a)(l-P)(u 2 =F V^«i) 

Reduce in ( |4.7] , |4.8[ ) hyperelliptic integrals to elliptic ones according 
to (pjpp. 



dx 



>/(l-x=«)(l-4 



+ 



dx 



a/(1 - x 2 )(l - k±x 



2 „2^ 



1t±, 



One can further express the symmetric functions of fix, fi 2 , Vi, v 2 on 
V x V in term of elliptic functions of tori E±. To this end we introduce 
the Darboux coordinates (see | 23fl , p. 105 ) 

Xi = sn(w + , /c + )sn(w_, £;_), 
(4.10) X 2 = cn(w + , k + )cn(u-, k-), 

X 3 = dn(u + , fc + )dn(w_, &;_), 

where sn(w ± , k±), cn(M ± , k±), dn (w±, A;±) are standard Jacobi elliptic 
functions. 

We apply further the addition theorem for Jacobi elliptic functions, 
sn(«i + u 2 , k) 
cn {u\ + u 2 , k) 
dn(ui + u 2 , k) 

where we denoted Sj = sn(«j, k), q = cn(«j, k),di = dn(«j, k), i = 1, 2 





- s 2 


s\c 2 d 2 


— S 2 Cidi 


s\Cid 2 


— s 2 c 2 <ii 


s±c 2 d 2 


— S 2 CiC?i 


Sid\C 2 


— s 2 d 2 S\ 


s\c 2 d 2 


- s 2 c\d\ 



and formulae (|3.7| , |3.9| ) for the Kleinian hyperelliptic functions. 
The straightforward calculations lead to the formulae 

(I - a)(l - {3)(a{3 + p 12 ) 



(4.11) 



x 2 = 



X, 



(a + P)(pi2 - a(3) + a(3p 22 + p n ' 
(1 + a/3) (a/3 - p 12 ) - a/3p 22 - p u 

(a + p)(pi 2 - aj3) + app 22 + p u 

af3p 22 - pn 

(a + P)(pi2 - a(3) + a(3p 22 + pn 
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The formulae ( 4.11]) can be inverted as follows 

x Mx 2 + x 3 ) - B(X 3 
(4.12) pn = ( B -l)^ ~ 



(4.13) pu = (B-l 

(4.14) p 22 



Xt + X 2 -1 

1 + X 1 -X 2 



x 1 + x 2 -r 

A(X 2 -X 3 )+B(X 3 -1) 



X 1 +X 2 -l 

where A = a + /3, B = 1 + a(3. 

The obtained results permit to present few solutions in elliptic func- 
tions of the initial problem, which are quasi-periodic in (. Using ( |4.13|) 
and ( f4.14|) for solutions of the (|2.1|) in the form ( [3.23| ) we have 



2 



-B-l 



1 f 2 A(X 2 -X 3 ) + B(X 3 -V 
1 + Xi-X 2 



1 at - a 2 V X x + X 2 -1 



X 1 + X 2 -l 



2 9 1 ( 2 A(X 2 -X 3 ) + B(X 3 -1) 

q2 ~ V 2 x l + x 2 -i a2 



[B-l 



1 + Xx - X. 



X t + X 2 -1 



where 

(4.15) u± = -2y/(l - a)(l - 0)(x =f c) 

and c is the constant depending on initial conditions. 

We also remark, that the derived quasi perodic solution was associ- 
ated with the Jacobi reduction case in which the ultraelliptic integrals 
were reduced to elliptic ones by the aid of second order substitution. 
This means on the language of two-dimensional ^-functions, that the as- 
sociated period matrix is equvalent to the matrix with the off-diagonal 
element T\ 2 = |. Such the reduction case was considered in various 
places (see e.g.flTj). Solutions of this type for nonlinear Schrodinger 
equation (a = 0) are recently obtained in HTT| . 

The anologous technique can be carried out for other well docu- 
mented case of reduction , when t\ 2 = 1/N and the N — 3,4, ... . In 
general such the reduction can be caried out for covers of arbitrary de- 
gree within the Weierstrass-Poincare reduction theory (see e.g. |29|, [7|). 

5. Elliptic periodic solutions 



In this section we develop a method (see also p7|, IS, 14|) which 



allows us to construct periodic solutions of fl2TT|) in a straightforward 
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way based on the application of spectral theory for the Schrddinger 
equation with elliptic potentials 0, |32|. We start with the formula 
(|3.14j ) and with equation with equation for Baker function \I/(A; u). 

dx 2 

where we identify the potential 



(5.1) — *(A, u) - U*(ar, u) = (A + -±)*(A, u), 



U = 2p 22 + 7"4- 


We assume, without loosing generality, that the associated curve has 
the property «4 = 0. To make this assumption applicable to the initial 
curve of the system fl2.1|) being derived from the Lax representation, 
we undertake the shift of the spectral parameter, 

2 2 

(5.2) A — >A + A, A=-ai + -o a . 

5 5 

Suppose, that U be two gap Lame or two gap Treibich- Verdier po- 
tential, what means, that 

N 

(5.3) U(x) =2^2p{x- Xi ), 

i=i 

where p(x) is standard Weierstrass elliptic functions with periods 2u, 2u' 
and numbers Xi takes the values from the set {0, u)\ = uj,u 2 = oj + 
= u'}. It is known, that the set of such the potentials is ex- 
hausted by six potentials^], [H]] 

= 6p(ar), 

= 6p(x) + 2p(x + uJi), i = 1,2,3, 

= 6p(x) + 2p(x + UJi) + 2p(x + Uj) , 
i^j = 1,2,3, 

- Qp(x) + Qp(x + tUi), i = 1,2,3, 
3 

= 6p(x) + 2^p(x + uj^, 

i=l 
3 

U 12 (x) = 6p(x) + 6^p(x + uj,i), 

i=l 

where the subscript shows the number of 2p functions involved and 
display the degree of the cover of the associated genus two curve over 
elliptic curve. Because the last three potentials can be obtained from 
the first three by Gauss transform we shall call the first three as basis 
potentials. The potential (|5.4|) is two gap Lame potential, which is 



(5.4) 


U 3 


:*) 


(5.5) 


U 4 


[x) 




U 5 


[x) 


(5.6) 








u 6 


[x) 




U 8 


[x) 
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associated with three sheeted cover of elliptic curve; the potentials 
(|5.19| . |5\6| ) are Treibich-Verdier potentials |46|, associated with four 
and five sheeted cover correspondingly. 

To display the class of periodic solutions of system ( j2.1|) we introduce 
the generalized Hermite polynomial T{x, A) by the formula 

(5.7) T{x, A) = A 2 - n 22 (x)X - n 12 (x) 

with n 22 (x) and iri 2 (x) given as follows 



N ^5 



3=1 



ni 2 (x) = -3 ^ p(x - Xj)p(x - Xj) - 



Ng 2 



i<j 



(5.8) 



i<3 \j=l 



where Xi are half-periods and N is the degree of the cover (see for 



example (JT8]] ) . The introduction of this formula is based on the possi- 
bility to compute the symmetric function fi\fi 2 in terms of differential 



polynomial of the first one with the help of the equation (|3.14|) , which 
serves in this context as the "trace formula" f49|j . 
The solutions of the system (|2.1|) are then given as 



(5.9) 



ql{x) 



2^' ai - A) , ql(x) 



JF(x, a 2 - A) 



The final formula for the solutions of the system fll.lp then reads 




l n J r (x } a 1 


-A) 


a\ — 


a 2 


n J r (x,a 2 


-A) 


a 2 - 






exp ^ %a\t — -v(a,i — A) 



dx 



i r 

exp < ia 2 t - -i>(a 2 - A) / 



^(x, a\ — A) 



da; 

T{x, a 2 - A) 



where we used ( [OP and ( |2.6|) . 

It is important to remark for our consideration, that if the poten- 
tial is known, then the associated algebraic curve of genus two can be 



described with the help of the Novikov equation [35]. Let us consider 
the two-gap potential for normalized by its expansion near the singular 
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point as 



6 

x 2 

Then the algebraic curve associated with this potential has the form 



(5.11) U(x) = — + ax 2 + bx 4 + cx 6 + dx 8 + 0(x 10 ) 



Fi • 7 I 2 ■ 7 

v 2 = \ 5 -^a\ 3 + —b\ 2 
2 2 

'3 4 • 7 o 3 3 • 11 \ . 3 4 • 17 . 3 2 - 11 • 13 



(5.12) + ^-—a 2 + ^^cj \-^^ab+ d. 

We shall consider below examples of genus two curves, which are 
associated with the two gap elliptic potentials (|5.4j ), ( |5.19|) and (|5.6|). 



Consider the potential U3 and construct the associated curve (|5.12 



(5.13) C 2 = (A 2 - 3# 2 )(A + 3ei)(A + 3e 2 )(A + 3e 



3/ 



The Hermite polynomiaal J r 3(p(x), A) |47|] associated to the Lame po- 
tential ( |5.4j ), which is already normalized as in ( |5.11| ) has the form 



(5.14) Fs(p(x), A) = A 2 - 3p(x)X + 9p 2 (x) - ~<? 2 . 



Then the finite and real solution to the system ( |2.1| ) is given by the 
formula ( [5.9D with the Hermite polynomial depending in the argument 
x + oj' (the shift in uj' provides the holomorphity of the solution). The 
solution is real under the choice of the arbitrary constants a 12 in such 
way, that the constants ai i2 — A lie in different lacunae. According to 
(|2.6| ) the constants Cj are then given as 



a 



o 4z/ 2 (a 4 -A) 



where A is the shift ( |5.2D and v is the coordinate of the curve ( |5.13|) 
and the integrals H and F have the following form 

25 ^ a2 ' T^ 2 ' 
F = ^ (ai + a 2 ) 3 - hjl - hjl - ^-g 3 - ^g 2 (a : + a 2 ) . 

These results are in complete agreement with solutions obtained in 
37f by introducing the ansatz of the form 



^Ap(x) 2 + Bip(x) + d, i = l,2 



with the constants A^Bi, Ci which are defined from the compatibility 
condition of the ansatz with the equations of motion. In the foregoing 
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examples we are considering the solutions of the form 

ft 0*0 = V Qi(p(x)), 

where Q, are rational functions of p(x). 

Consider with this purpose the Treibich Verdier potential 

(5.15) \J A (x) = 6p(x) + 2p(x + u 1 )-2e u 

associated with four sheeted cover. The potential is normalized accord- 



ing to (|5.11|) . The associated spectral curve is of the form 



(5.16) v 2 = 4(A + 6 ei ) J](A-A fc ) 



k=l 



(5.17) A 1)2 = e 3 + 2e 2 ±2y/(5e 3 + 7e 2 )(2e 3 + e 2 ) 
A 3 ,4 = e 2 + 2e 3 ± 2^(5e 2 + 7e 3 )(2e 2 + e 3 ), 

The Hermite polynomials are given by the formula 

(5.18) F(x, A) = A 2 - (3p(x) + p(x + uji) - e x )A 

+ 9p(x)(p(x) + p(x + u) — ei) - 3eip(x + Wi) 
9 2 



The finite real solution of ( |2.1j ) results the substitution this Hermite 
polynomial T{x + lu', A) into ( |5.9| ) depending in shifted by imaginary 
half period argument into the formula (answer). To provide the reality 
of the solution we shall fix the parameters cij — A in the permitted 
zones. The constants Cj are computed by the formula (j|) at which v 
means the coordinate of the curve (|5.16j ). 

Consider further the Treibich Verdier potential 

(5.19) U 5 (x) = 6p(x) + 2p{x + uo 2 ) + 2p{x + cu 3 ) + 2e u 

associated with four sheeted cover. The pontial is normalized according 
to ( |5.11| ). The associated spectral curve is of the form 

v 1 = (A + 6e 2 -3e 3 )(A + 6e 3 - 3e 2 ) 

x [A 3 + 3eiA 2 - (29e 2 - 22e 2 e 3 + 29e 2 )A 



(5.20) + 159(e 3 + e 3 ) - 51e 2 e 3 (e 2 + e 3 ) 

The associated Hermite polynomials are given by the formula 

F(x,\) = A 2 - (3p(x) + p(x + u 2 )) + p(x + u 3 ) + ei)A 

+ 9p(x)(p(x) + p(x + u 2 ) + p(x + u 3 )) + 3p(x + u 2 )p(x + uj 3 ) 

39 

+ 3ex(3p(x) + p(x + lo 2 )) + p{x + uj 3 )) - —g 2 + 54e 2 . 
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The solution of the system results the substitution of these expressions 
to ( |5.9|) as before, but this solution is blowing up. 

We remark, that since McKean, Moser and Airlault paper is well 
known, that all elliptic potentials of the Schrodinger equations and 
their isospectral transformation under the action of the KdV flow has 
the form, 

N 

(5.21) U(x) = 2^p(x-x J (t)), 

i=i 

The number N is a positive integer N > 2 (the number of "particles" ) 
and the numbers x — (xi(t), . . . , x^(t)) belongs to the locus C^, i-e., 
the geometrical position of the points given by the equations 

(5.22) C N =\(x);J2 p'Mt) ~ X M = 0, j = 1, . . . iV 

If the evolution of the particles Xi over the locus is given by the equa- 
tions, 

then the potential ( |5.21| ) is the elliptic solution to the KdV equation. 
Henseforth the elliptic potentials discuss can serve as imput for the 
isospectral deformation along the locus. Moreover these elliptic poten- 
tial do not exhausted all the viriety of elliptic potential; we can mention 
here the elliptic potentials of Smirnov (|40], ^T|) for which the shifts Xi 
are not half periods. The involving of these objects to the subject can 
enlarge the classes of elliptic solutions to the system fll.lp 



6. Conclusions 

In this paper we have described a family of elliptic solutions for the 
coupled nonlinear Schrdinger equations using Lax pair method and the 
general method of reduction of Abelian functions to elliptic functions. 
Our approach is systematic in the sense that special solutions (periodic, 
soliton etc.) are obtained in a unified way. We considered only the 
family of elliptic solutions associated with the integrable case 1:2:1 
of quartic potential, the approach developed can be applied to other 
integrable cases being enumerated in the introduction. 

In fiber optics applications of the periodic and quasi-periodic waves 
are of interest in optical transmission systems. 
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